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Abstract 



(^ , The system of four point vortices in the plane has relative equilibria 

I*^ ' that behave as composite particles, in the case where three of the vor- 

tices have strength — r/3 and one of the vortices has strength F. These 
0^ . relative equilibria occur at nongeneric momenta. The reduction of this 

^N ' system, at those momenta, by continuous and then discrete symmetries, 

(^ I classifies the 4-vortex states which have been observed as products of col- 

O |. lisions of two such composite particles. In this article I explicitly calculate 

these reductions, and show they are qualitatively identical one degree of 
freedom systems on a cylinder. The flows on these reduced systems all 
C^ ' have one stable equilibrium and one unstable equilibrium, and all the or- 

bits are periodic except for two homoclinic connections to the unstable 
equilibrium. 

In the system of TV vortices in the plane the n}^ vortex has location 
Xn + iyn € P^^ = (C^)^. The Hamiltonian and symplectic form are 



Zrn 



1 ^ 

i?^'--^^r„r,„in|z„-z„|2, c^pi = 0r„^o, (i) 



n=l 



where ti;o(a, 6) = — Im(a6) for complex numbers a, 6 G C. This system admits 
the symmetry group SE{2) — {(e'^,a)} of Euchdean symmetries acting diago- 
nally on each factor C of P^^ by (e*^, a) ■ z = e^^z + a. A momentum mapping 
is 
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where se(2)* is identified with se(2) = {{0,a)} = R'^ by the standard inner 
product of R"^. For more details on the system of N vortices on the sphere or 



plane, see Kidambi and Newton |1998], Pekarsky and Marsden | 1998( |, and the 
references therein. 

In Patrick 1 2000( | I have shown that the point of phase space corresponding 
to a central vortex of strength F at the origin surrounded at distance a by 
three symmetrically placed outer vortices of strength — r/3 is a formally stable 
relative equilibrium for this system; that is 
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is a formally stable relative equilibrium of the planar point vortex system with 
TV = 4 and 
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The momentum mapping is equivariant for these vortex strengths because the 
sum of the vortex strengths is zero. The relative equilibrium itself does not 
translate but merely rotates with generator Oe = T/Sna^, and by direct substi- 
tution the momentum is ^e = (rQ!^/2, 0). This momentum value is nongeneric, 
since it has isotropy group all of SE{2), and the position of this relative equilib- 
ri um is un s table under small perturbation to nonzero momentum: as discussed 
in Patrick [^000 1, the relative equilibria will move about on the plane as a com- 
posite particle. To date these are the only relative equilibria of point vortices 
which are both formally stable and have nongeneric momentum. The particular 
system of four vortices with the strengths F^ is singled out by this fact. 

When two such relative equilibria collide the result may be localized 4- vortex 
states far from such relative equilibria but each near to momentum /ig, and these 
too have been observed to move about on the plane as composite particles. The 
phase space of classifying these states, since they are far from the relative equilib- 
rium (H), yet still close to momentum /ig, will be the reduced phase space of the 
4- vortex system at momentum /ig . This phase space has dimension 4x2 — 6 = 2 . 
Here I calculate this phase space and, qualitatively, its associated Hamiltonian 
flow. One can find other reductions on various point vortex systems, in the 
sphere and in the plane, in Adams and Ratiij | 1988|, Aref and Pomphrcy |1982|, 



Eckhardt and AreJ 



and Pekarsky and Marsden [1998|. The general is 



sue of drifting relative equilibria at nongeneric momenta is discussed at length 
in the sequence of articles Patrick [1992|, Patrick [1996|, Patrick |199E], and 



Patricli 11999] 



To calculate the quotient 

one may first translate the central vortex to the origin, whereupon, by (^, 
the sum of the outer vortex positions must vanish. The resulting 4-dimensional 
vector space C^ — {(wi +ivi, U2 + iv2)} is spanned by the following orthonormal 



basis {Ei} of outer vortex positions: 
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or more specifically the column vectors of the matrix 



Zl 

A 



E4 



' Z2 



E 



V3 



1 





1 








1 





1 


1 

2 


2 


1 

2 


2 


2 


1 

2 


2 


1 

2 


1 
2 


71 
2 


1 
2 


_V3 
2 


V3 
2 


1 

2 


2 


1 

2 



The SE{2) symmetry then remaining on C^ is 50(2) acting by diagonal clock- 
wise rotation, and since the above basis is orthonormal, the level set of Jp' 
becomes 

{ (mi + ivi,U2 + iv2) ■■ ui^ + wi^ + -^2^ + W2^ = 3a^ }. 

The Hopf variables 

2 —1 

Wl = -;-^{uiV2 - U2V1), W^ = ^-^((Ul^ + Wl^) - ("2^ + 1'2^)), 

2 1 

^2 = -—^{uiU2 +1^11^2), W4 = ^-^(("1^ +Wl^) + ("2^ +""2^)), 

with image the semi-algebraic set 

Wl^ + W2^ + ^3^ — ^^4^ = 0, W4 > 0. 

are a quotient map for the diagonal 50(2) action on this space, so restricting 
these to the w^ — 1 sphere gives the quotient as 

(^"^e ={{W1,W2,W3) ■.Wi^+W2'^+W3^ = 1}. 

One partial section to the Hopf mapping, obtained by imposing vi — 0, is given 

by 
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and a section to the map from PP' to C^ may be obtained by setting Z4 — and 
then multiplying the vector (ui, wi, M2, V2) by E to obtain the triple (zi, Z2, ^3)- 
Using this section to pull back the symplectic form on PP' gives the reduced 
symplectic form 
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UJS2{w){wi,W2) = -w ■ {Wi X W2). 



To summarize, the symplectic reduced phase space is the 2-sphere with sym- 
plectic form the constant multiple Ta^/A of the standard symplectic form on 
the 2-sphere obtained by left-hand reduction of the cotangent bundle of 50(3). 
The action of the discrete group 63 corresponding to permutation of the 
identical outer vortices descends to the reduced space since it commutes with 
the action of SE{2). Using the above section this reduced discrete action is 
easily calculated to be the restriction to the 2-sphere of the linear representation 
cr : S3 -^ 50(3) given by 
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The squares of the interparticle distances can be descended to the reduced 
space, giving a^/2 times each of the 6 linear functionals 

Z12 = -3(\/3wi - W2 - 2), ^13 = 3(V3wi + ?«2 + 2), ^23 = -6(^2 - 1) 

Z4I = 2{W2 + 1), ?42 = -(V3wi +11^2- 2), ^43 = \/3wi - W2 + 2 

and in terms of these the reduced Hamiltonian is 

a^\ (^1^42^43)^1 




^12^13^23 



The linear functionals Uj are independent of ws, and are zero for the two equi- 
lateral triangles circumscribing the unit circle in the {wi, W2) plane (see the left 
of Figure (0)). The Hamiltonian is finite except for the points in this plane 
where the triangles touch the circle, where it is positive infinity on the triangle 
that points downward, corresponding to collision states of two outer vortices, 
and negative infinity on the triangle pointing upward, corresponding to colli- 
sions states of an outer vortex with the central vortex. The Hamiltonian has an 
amusing relationship with the geometry of the collision triangles: it differs by 
a constant from the function that is the logarithm of the product of the three 
distances to the — 00-collision triangle, cubed, divided by the product of the 
three distances to the +cx3-collision triangle. 
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H = 00 



Figure 1: Left: the level sets of the linear functionals lij are the triangles shown. Right: The 
reduced Hamiltonian i/^^ ; the sphere has been stereographically projected from the south pole 
and then the plane compressed to a square, so that the boundary of the square is the south pole. 
Clearly visible are the ±00 collision states and at the center the stable equilibrium. Another 
stable equilibrium was, on the original sphere, located at the south pole, so it occupies the 
boundary of the square. Between the +00 collision "mountains" and the stable equilibrium 
at the center are three saddle points and there are three other saddle points between the +00 
collision "mountains" and the edge of the square. 

The equilibria, and hence the relative equilibria of the unreduced system, 
may be found by calculating the critical points of iJ^^ . The south pole is the 
stable equilibrium corresponding to the formally stable relative equilibrium (^ , 
and the north pole corresponds to the mirror image of this, obtained by ex- 
changing Z2 and Z3 in (j^). There are six unstable equilibria, namely the points 
of the orbit of the equilibrium wi = 0, W2 = v3 — 1 under the action of 63. 
One relative equilibrium corresponding to these unstable equilibria is 

zi = a\^, z2^- — \^'] p(\/3-3)i, 23 = ^2, 2:4 = 0. (4) 

^ 2v2 

A graph of the reduced Hamiltonian as a function of wi and 1/72 is shown on the 
right of Figure (|l]). 

The above is a complete analysis of the reduced space and the phase portrait. 
However, to find the most succinct classification of the localized 4- vortex states, 
one should also quotient by the action of 63. Then the unstable equilibria will 
collapse to a single point, the two stable equilibria collapse to a single point, 
and there will be one +00 collision state and one —00 collision state. It is easy 
to calculate that the stable equilibria have isotropy the cyclic group 33, the 
collisions have isotropy the cyclic group 32, and these are the only points of the 
reduced phase space with isotropy. 

To find the quotient of the reduced space by 63 one can find its image under 
sufficiently many functions invariant under the action of 63. By averaging up 
to degree 4 polynomials in wi, W2,W3, one finds the invariant polynomials 

Pl = U)3^, p2 = W2{3wl-wl), 

P3 = WiW3{3w2^ - Wl^), Pi = Wl^ + W2^ + Ws"^ 

and the relations 

Pi{{pi - Pif + P2^) + P3^ = 0, 0<pi<p4. (5) 





Figure 2: Left: top half of the semi- algebraic set m). The +00 collision state is to the front 
and the —00 collision state is behind and obscured; the singularity at (pijPiJJs) = (1, 0,0) is 
the relative equilibrium M). Right: same as left but deformed according to ([Tl). The resulting 
surface may be radially projected to the sphere and then stereographically projected to the 
plane from the point (51, 92,53) = (0, —1,0). 



The map (wi,W2,W3) 1-^ (pi,p2,P3,P4) can be seen as the composition of simpler 
steps, namely 



Pi = W3, 
followed by 



P2 + ip3 = {—W2 + iwi)^ 



pi = Wl^ 



W2^ + W3^, 



Pi =P1 



P2=P2, P3=PlP3, Pl=Pi, 



and hence it is easily seen that (wi,W2,W3) 1-^ {pi,p2,P3,P4), with image the 
semi-algebraic set (||), is a quotient map for the action of 63 on C^. Thus (||) 
restricts to a quotient map of the reduced space to the semi-algebraic set 



Pl{{pi - If + P2^) + P3^' 

The Hamiltonian this reduced space is 



0, 0<pi<l. 
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The set (O), the top half of which is graphed on the left of Figure (0), is 
homeomorphic to a two sphere and has three singularities, corresponding to 
the two collision types of the reduced phase space and to the stable relative 
equilibrium (H). If the ±oo-collision states are removed then the phase space 
is a cylinder; however, the symplectic form of this phase space is undefined at 
one point, namely the stable equilibrium corresponding to the relative equilib- 
rium (0) . An explicit map to this cylinder can be constructed by deforming (m 
according to 



91 = Pi - ^(1-P2^), 



92=P2, 93=^3, 



(7) 



resulting in the surface plotted on the right of Figure (0) . The collision states 



(Pl,P2,P3) = (0,±l,0) 





Figure 3: Right: the result of projecting, by M), simulations of the ^-vortex system to the 
toral reduced space as described. The horizontal variable is d^^ and the two vertical sides 
of the phase portrait should be identified. Left top to bottom: motions of the vortices on 
the plane near the +00 collision, the top half of the homoclinic orbit, the bottom half of the 
homoclinic orbit, and the —00 collision state. The pictures on the left are reconstructions of 
orbits on the right. Left, the third graph from the top: the motion of the vortices is a close 
pair for strengths —1 orbiting two vortices, one of strength 3 and the other of strength —1, 
and periodically the vortex on the inside is exchanged with one of the vortices on the outside; 
this corresponds to the lower half of the homoclinic orbit on the right. Left, second graph from 
the top: the same configuration but the two outer vortices making up the close pair exchange 
positions; this corresponds to the upper half of the homoclinic orbit on the right. 
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Figure 4- The space of localized 4-vortex states is the "Y" shaped union of the two lines 
above. The stable relative equilibrium (M) is represented by the endpoint of the line branching 
upwards towards the right. At the top and bottom endpoints of the vertical branch are the 
±00 collisioristates respectively. The bifurcation point in the middle is the unstable relative 
equilibrium (Uj. From that point proceeding upwards immediately yields to modes of motion 
corresponding to the two connecting homoclinic orbits, as indicated on the left. 



are unchanged but the result of the deformation uniquely intersects radial 
half lines from the origin, and thus can be radially homeomorphed to the 
sphere. Then, stereographically projecting from the coUision state (91, 92, 93) = 
(0, —1,0) results in the plane C \ {0}, and the pairing of the angle 0^^ at the 
origin together with the logarithm /i^^ of the radius results in a map to the 2- 
torus such that the ±00 energy collision states are mapped to ±00 respectively. 
The final map to the cylindrical phase space {{h, 9)} is then 

" = arctan I — I , (8) 




and Figure (||) shows the resulting phase portrait. 

As stated at the beginning, the prime motivation of this article has been 
the classification of 4-vortex states, produced through collisions of two stable 
relative equilibria (3), with momenta near /ig = (ra^/2,0) but far from the 
relative equilibria (3). Such localized states ought to be roughly regarded as 
orbits of the phase space in Figure (H), since there will be a fast time scale 
roughly corresponding to motion along those orbits and the localized states will 
present external properties corresponding to averages over that fast time scale. 
Thus, the space of localized 4-vortex states corresponds to the quotient space 
of the torus by the orbit equivalence relation of the phase space of Figure (0) . 
This quotient is shown in Figure (0). 
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